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Quenching effects in strong penguin contributions to e'/ e 

Jack Laiho" 

"Department of Piiysics, Princeton University, Princeton, NJ 08544 

Quenching effects in strong penguin matrix elements are investigated. A lattice determination of a^^, 
constant that appears in the quenched ChPT relevant for the lattice analysis of if — > tttt matrix elements, shows 
that this constant is large. The original RBC analysis of Qa matrix elements is revisited in light of this result. 
Also, the numerical effects of choosing the singlet Golterman-Pallante method of quenching Qe is investigated. 



the 



> 



0^ 
O 

O 

CD 



X 



1. Introduction 

There have been several recent lattice attempts 
to calculate Re(e'/e), the direct CP violating pa- 
rameter in K TTTT decays. These include the 
attempts using domain wall fermions by the CP- 
PACS [T] and RBC 2^ collaborations. A notable 
feature of both of these calculations is that their 
central values differ drastically from experiment, 
though the approximations made were rather se- 
vere. In particular, the quenched approximation 
and leading order chiral perturbation theory were 
employed in both calculations. 

Although the quenched approximation is un- 
controlled, where quenched lattice results have 
been compared to experiment, at least for sim- 
ple quantities such as masses of flavored mesons 
and decay constants, the agreement is at or better 
than ~ 25%. However, there is no good reason 
to think that this agreement should hold for all 
low-energy hadronic phenomena. In fact there is 
a particular difficulty with quenched e'/e that is 
not present in other quenched calculations. 

Since the original RBC e'/e result was reported, 
it was shown by J3 _4_ that there is an ambiguity in 
defining the quenched strong penguin four-quark 
operators. This ambiguity affects the operator. 



Qe = s„7^(l - 7^)4 ^ g67''(l + 7^)'7a. 



(1) 



whose K — > TTTT matrix element has a large con- 
tribution to e'/e. There are at least two ways of 
defining this operator in the quenched theory, and 
we dub the calculations done with these different 



definitions the quenched singlet and non-singlet 
methods. The non-singlet method corresponds to 
the choice of the operator made in the standard 
(original RBC) approach, while the singlet choice 
corresponds to the Golterman-Pallante proposal. 
It is not known a priori whether either of these 
two methods gives results close to those of the 
full theory, and debates on this matter are likely 
to be resolved only when we have the results of 
dynamical calculations with the physical number 
of light flavors. It is possible that an appropriate 
linear combination of the two methods may give 
a reasonable approximation to the full theory; see 
the contribution of Norman Christ to these pro- 
ceedings for such a proposal. 

For the non-singlet method used (implicitly) in 
RBC's previous work, it was shown that K — + tttt 
may have a large contribution from an additional 
low energy constant (LEC), as suggested by the 
large Nc approximation, where Nc is the number 
of colors [Sj. Since for all AI — 1/2 amplitudes, 
the leading order ChPT analysis employed in |2] 
was that of the full theory, it seems desirable to 
revisit the analysis of Qe matrix elements in light 
of the new results in quenched ChPT. We em- 
phasize that this ambiguity is present only for 
the strong penguin operators, and so does not af- 
fect the Qs contribution to e'/e. This quenching 
ambiguity also does not affect Re Aq, Re A2 and, 
therefore, the AI =1/2 rule. 

2. Singlet and non-singlet methods 

We discuss the situation for Qg in order to il- 
lustrate the subtlety involved in quenching the 
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strong penguin operators. From Eq (1) it can be 
seen that the right part of Qg is a sum over ac- 
tive Ught flavors, so that in the full theory the 
right part is a flavor singlet under the symmetry 
group SU(3)fl. In the quenched theory one has at 
least two options. In the first option one chooses 
to sum over all the quarks, including valence and 
ghost. In this case the right component of the op- 
erator transforms as a singlet under the extended 
symmetry group; this is called the quenched sin- 
glet option. 

In the second option, one may choose to sum in 
Eq over only the valence quarks. In this case 
the operator can be decomposed into two opera- 
tors, one of which transforms as a singlet under 
the irrep of the extended symmetry group, while 
the other does not (rather, for Qe, it transforms 
in the adjoint representation); we choose to call 
this the quenched non-singlet method. In this 
method a new constant, , appears at lead- 
ing order (LO), and a large Nc estimate |S] has 
suggested that this constant may be large. 

3. Measuring a^^ 

It is difficult to numerically disentangle the 
a^^ log term from the linear term and possible 
higher order effects where it appears in K ^ 
It was for this reason that 0] provided a ma- 
trix element to which a^^ contributes at 0(p*'), 
so that it may be obtained more readily in the 
chiral limit. This can be done if one considers a 
matrix element where the ghost quarks can ap- 
pear on external lines. 

Golterman and Pallante proposed looking at 
the matrix element, K ^ 0, where a^^ can be 

obtained in the chiral limit. The A" is a me- 
son made of a quark anti-ghost. We have shown 
that there is a graded CPS symmetry that is rel- 
evant for matrix elements with ghosts in external 
states that is analogous to the usual CPS sym- 
metry used to constrain the power divergences of 
kaon decay amplitudes. Using this symmetry, we 
have shown that the K matrix element has 
a quadratic divergence. Also, because the chi- 
ral symmetry of domain wall fermions with finite 
fifth dimension is not exact, there is a contamina- 
tion in the chiral limit of 0{mres/o,^) where rrires 




Figure 1. mQ^'''\K) + 

{0\Q^^^\K))/{0\sj5d\K) as a function of 
irid (big circles) and nid + rrires (squares), and 
/V(2a(3.3))(^|QQwS|^^ (small circles) as a 
function of rUd equal to the degenerate quark 
mass TO/. The pre- factor of {n\Q^^'^\K) was 
chosen so that the two methods would have the 
same value in the chiral limit according to ChPT. 

is the contribution to the fermion mass from chi- 
ral symmetry breaking. Figure 1 shows that the 
size of this ambiguity is of the same order of the 
chiral limit we are trying to extract. 

For this reason we have proposed an alterna- 
tive method for obtaining that does not suffer 
from this problem. Using the graded CPS sym- 
metry we have shown that for degenerate quark 
masses, the amplitude isT — > tt does not have di- 
vergences. This matrix element calculated on the 
same configurations is plotted in Figure 1, and 
the value of a^^ so obtained is in rough agree- 
ment with the large Nc calculation implying 
that this constant must be taken into account in 
the RBC analysis. 

4. EfTects of a^^ in the non-singlet analysis 

The LO ChPT prediction for the K matrix 
element including the a^^ does not produce a 
good fit to the data (uncorrelated x^/dof « 30), 
and it is possible that higher order effects need 
to be included in the analysis. The data appears 
quite linear, and cannot accommodate the large 
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log in the analytical expression. Unfortunately, 
the NLO contributions for the non-singlet method 
are not currently known, so it is not possible to do 
a systematic analysis. The corrections to K ^ 
for the singlet part of the Qe operator required 
a one- loop calculation to go to NLO [0{p'^)] [7]. 
The non-singlet contribution to 0{p'^), however, 
would require a two-loop calculation, but this is 
not likely to be available anytime soon. 

It is clearly difficult to separately resolve non- 
linearities using the numerical data if there is, 
in fact, a cancellation between logarithmic and 
quadratic terms over the range of quark masses 
used in our simulations. Although there is reason- 
able agreement between the results for the sub- 
tracted K TT matrix element leaving out the 
heaviest mass data points in {0\Qe\K'^) versus in- 
cluding a quadratic term in the fit and using all 
ten data points, we do not believe we have enough 
control over the subtraction in order to obtain the 
needed LEC's from K ^ n. We, therefore, can- 
not claim a result for K vrvr in the non-singlet 
method. 
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Figure 2. The lower data set is the matrix element 
{n~^\QQ^\K^)sub where the divergent piece has 
been subtracted. This was obtained using the 
singlet method, and the fit is to the known NLO 
form. The upper data set is the original RBC 
result for {TT'^\Qe\K~^)sub with a linear fit shown 
for comparison. 

NLO formulas as the result can be affected by a 
factor of two or more. 



5. Results for the singlet method 

Although the quenched uncertainties make it 
unclear as to which method is best to approxi- 
mate the full theory, the singlet method has the 
practical advantage that the NLO formulas are 
known, making it a useful test of the whole for- 
malism, including a comparison to (quenched) 
continuum methods. The quenched NLO formu- 
las do indeed produce good fits to the data. We 
are able to obtain a subtracted K —> n amplitude 
which is shown in Fig 2, along with the original 
RBC data reported for subtracted if — > tt (where 
the a^^ term was not taken into account and a 
linear fit was used), for comparison. The slope 
of this plot as a function of my gives the LEG 
needed for LO K —^ tttt determinations, and one 
can see that the slope is enhanced by roughly a 
factor of two in the singlet method over the orig- 
inal RBC result. On the other hand, a complete 
NLO fit shows that the linear part of the singlet 
fit (the slope in the chiral limit) is roughly par- 
allel with the slope quoted for the original RBC 
data. Thus, it is extremely important to have the 



6. Outlook 

The large quenching effects show that dynam- 
ical calculations will be necessary in order to sig- 
nificantly reduce this large systematic uncertainty 
in e'/e. Also, if the reduction method of [HI is 
used, it appears necessary to have NLO ChPT in 
order to extract the needed LEC's from lattice 
data. 
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